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ABSTRACT 



The point-spread function (PSF) of an adaptive optics (AO) system is often poorly known. 
2. . This ignorance can lead to significant systematic errors. Since the degree of AO correction is 

in ' sensitive to the observing conditions: seeing, wind speed, brightness of the wavefront reference, 

etc., it would be desirable to estimate the PSF from the data themselves rather than from 
observations of a PSF star at another time. We have developed a method to estimate the PSF 
delivered by an AO system in the case where the scene consists of a crowded star field. We 
model the modulation transfer function (MTF) of several key components of the imaging system 
(atmosphere filtered by an AO system, telescope pupil, and pixel array). The power spectrum of 
the image, even a dense star field, can be used to constrain our model, which in turn can be used 
to reconstruct the PSF. In the case of circularly symmetric PSFs, we demonstrate that the power 
spectrum of the source distribution function can be successfully removed from the measured MTF 
and that our fit successfully recovers input parameters from a model data-set constructed from 
these parameters. We also show that the method yields reasonable fit parameters and a useful 
approximation to the PSF when applied to data from the laser guide star (LGS) AO system at the 
Keck Observatory. Comparison of Keck LGS/AO data and Hubble Space Telescope observations 
with NICMOS show that photometric accuracy of a few percent can be achieved for data with 
Strehl ratios as low as 4%. 

Subject headings: atmospheric effects — instrumentation: adaptive optics methods: data analysis — 
techniques: photometric — stars:supergiants: variable 

1. The Adaptive Optics PSF Wizinowich et al. 2004; Rousset et al. 2000). 

There is a strong incentive to develop and deploy 
Adaptive optics provides the ability to achieve A0 gystems becauge the cogt of a larg6i ground _ 

near diffraction-limited performance on large based telescope equipped with AO is modest com- 

ground-based telescopes (e.g., Herriot et al. 2000; pared fco that of & comparable aper ture space tele- 

T~ " 7 „ TO „ , scope. Most large observatories already have AO 

"'Now at the Department of Physics &c Astronomy, . 

UCLA, Los Angeles, CA 90095 systems either m place or under development, and 

1 Based on observations obtained at the W.M. Keck Ob- the data Ir0m these systems are beginning to yield 

servatory. scientific results (e.g., Carson et al. 2005; Minowa 



ct al. 2005; Melbourne et al. 2005; Max et al. 2005; 
Perrinet al. 2004). 

While AO images show clear advantages over 
seeing-limited data in resolution and sensitivity, 
there remain significant obstacles to the full anal- 
ysis of these data. Characterization of the AO 
point-spread function (PSF) is arguably the most 
serious of these challenges. Conventional AO sys- 
tems deliver modest correction (Strehl ratio, SR < 
0.5) over a field of view of a few tens of arc sec- 
onds (e.g., Fitzgerald & Graham 2006) — although 
some high-order natural guide star systems afford 
better performance (e.g., Carson et al. 2005). Ex- 
perience shows that the chance of finding a suit- 
able PSF star within this field is small. Moreover, 
because the precision and accuracy of wavefront 
correction is sensitive to the observing conditions, 
including seeing, wind speed and brightness of the 
wavefront reference it is often a futile exercise to 
observe a PSF calibration field at a different time. 

Changing atmospheric conditions and AO per- 
formance lead to variable correction, causing vari- 
ations in the amount of energy in the diffraction- 
limited core of the PSF relative to the uncorrected 
seeing halo. These variations affect the Strehl ra- 
tio. Systematic errors in photometry will occur if 
these Strehl variations are not taken into account 
in evaluating the photometric curve of growth and 
the corresponding aperture corrections. The angu- 
lar scales of diffraction and seeing are X/D te i and 
A/ro , respectively, where A is the observing wave- 
length, D te i is diameter of the telescope, and 7*o 
is the Fried parameter. For a large telescope at a 
good site, operating in the near infrared (IR), the 
ratio Dtei/ro ~ 10 — 20. Thus, the surface bright- 
ness of the seeing halo can be more than two or- 
ders of magnitude fainter than the PSF core, ren- 
dering it challenging to measure the seeing halo 
directly — even if a suitable star is present. The 
direct approach to estimating the PSF is further 
frustrated because science targets selected for AO 
observations will tend to have multiple sources per 
seeing disk. 

An innovative approach uses telemetry from the 
AO system to estimate the PSF from the resid- 
ual wavefront errors (Veran 1997). However, this 
method has not been implemented for all AO 
system architectures, and it cannot correct for 
non-common path wavefront errors between the 
wavefront sensor and the science camera. My- 



opic deconvolution is another promising alterna- 
tive, which simultaneously estimates the source 
distribution and the PSF (Jefferies & Christou 
1993). In the "StarFindcr" code (Diolaiti et al. 
2000), which is designed for analysis of variable 
Strehl-ratio observations of stellar fields, the PSF 
is extracted directly from the image frame. 

These approaches use information about the 
imaging system to various degrees. For example, 
in myopic deconvolution it is possible to specify 
the support of the PSF in frequency space (Jef- 
feries & Christou 1993). It is our objective to 
capture some rudimentary properties of the atmo- 
sphere, AO system, telescope and detector to build 
a model PSF, which can be constrained by AO 
observations of astronomical scenes. Our primary 
interest is in estimating the photometric curve of 
growth, so we consider only the azimuthally aver- 
aged PSF. Initial astronomical application of this 
algorithm was made by one of us (JRG) in de Pa- 
ter et al. (2004). 

In § 2 we present our description of an AO sys- 
tem, including partial wavefront correction (§ 2.1). 
In § 3 we describe how to interpret the power spec- 
trum of an AO image in terms of the MTF and 
the source distribution power spectrum. Section 4 
describes the influence function of the dcformable 
mirror in the Keck AO system. We use the re- 
sulting model to simulate Keck AO data (§ 5) and 
show that the PSF can be recovered from synthetic 
images (§ 5.1), and the results of application to 
Keck LGS AO data are outlined in § 5.3. 

2. The Transfer Function, Point Spread 
Function and Encircled Energy 

The complex optical transfer function (OTF) 
describes the response of a linear optical system 
at some spatial frequency in the image plane, v. 
The point-spread function (PSF) forms a Fourier 
pair with the system OTF. Each independent com- 
ponent of an optical system, from the atmosphere 
through to the detector, has its own OTF. The 
system OTF is the product of each of the sepa- 
rate OTFs (e.g., Schroeder 2000). For the bulk of 
this paper, we consider only circularly symmetric 
PSFs. In this case, the complex part of the OTF is 
zero and the system is fully described by the real 
part of the OTF, which is known as the modula- 
tion transfer function (MTF). We consider three 



principal components that contribute to the sys- 
tem MTF: the atmosphere (partially corrected by 
an AO system), the telescope pupil, and the de- 
tector. Our notation is summarized in Appendix 
A. 

2.1. Partial Wavefront Correction 

Suppose the phase in the pupil plane of the tele- 
scope at a point r is <j>(r) and that the phase cor- 
rection presented by the deformable mirror (DM) 
at this point can be described by the convolution 



<t>DM{*) 



[<j>(r')+S<j)(r')}h(r-r')dr', (1) 



where S<p(r) represents the error in the wavefront, 
and /i(r) is the influence function of the DM, i.e., 
the shape on the DM created by pushing a single 
actuator. The term <5</>(r) can be used to repre- 
sent a static wavefront aberration, or the wave- 
front measurement error. The corrected phase af- 
ter the DM is 



</>Ao(r) = <j>(r) - (f>DM(r). 



(2) 



Taking the Fourier transform of Eq. (2), substitut- 
ing <f>DM from Eq. (1), and applying the convolu- 
tion theorem, the expectation value of the power 
spectrum of residual phase errors at spatial fre- 
quency in the pupil plane, k = \k\, is 

|$ ao (k)| 2 = [1 - H{n)f |$(«)| 2 + H 2 ( K )Al , (3) 

where H (k) is the Fourier transform of the in- 
fluence function and A? is the variance of wave- 
front measurement error. In deriving Eq. (3), 
we assumed that S(f>(r) represents uncorrelated, 
zero-mean noise. The notation |<I>| 2 denotes the 
atmospheric phase error power spectrum (see 
§ A.l). The AO system can be thought of as high- 
pass filter in the pupil plane, with filter function 
[1 - H(k)} 2 (see § 4 and Fig. 1 (a)). 

The AO-corrected atmospheric structure func- 
tion, Dj, AO (r), can be found from Eq. (A4) by sub- 
stituting the filtered phase error power spectrum 
|$Ao(t)| 2 fromEq. (3). Since |$ao(«0| 2 and H(k) 
are both circularly symmetric the 1-d form of the 
integral can be used. The expression for the MTF 
of an AO-corrected atmosphere is then given by 
Eq. (A7), 



The final system MTF is given by the product of 
the MTFs for the atmosphere, the telescope (Eq. 
A8), and the detector (Eq. A9) 

T sys {v) = T0 AO (^) T pup (v) T plx (v). (5) 

The final step involves finding the PSF (Eq. A10) 
or the encircled energy (EE) (Eq. All) from the 
Hankcl transforms of T sys (y n ). 

3. Extracting the MTF from Astronomi- 
cal Data 

We now describe a procedure for estimating 
the MTF from observations of a crowded stellar 
field. Consider an astronomical scene consisting 
of a distribution of sources of amplitude and po- 
sition given by the function s(x). The resultant 
image, z(x), 



i(x) = s(x) * p(x) +n(x). 



(6) 



is the convolution of this scene with the PSF plus 
additive noise n(x). We assume that n(x) is zero- 
mean Gaussian noise. This assumption will re- 
quire our data are accurately sky-subtracted be- 
fore analysis. Taking the Fourier transform of 
Eq. (6) and applying the convolution theorem, we 
find that the power spectrum of the image, ^(f) | 2 , 
is the product of the source distribution power 
spectrum, \S(i>)\ 2 , and the PSF power spectrum, 
plus a constant noise floor, |iV| 2 . Since the PSF 
and the MTF are Fourier pairs, the power spec- 
trum of the image is given by 



\I(u)f = \S(v)\ 2 -\T(u) 



\N\ 



(7) 
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If the source distribution is spatially uncorrelated 
then IS^f)! 2 can be approximated as a multiplica- 
tive constant, and we should be able to estimate 
the MTF directly from |/(i/)| 2 . Using a model for 
the imaging performance of the AO system (e.g., 
§ 4), we may fit \I(u)\ 2 , yielding the MTF and 
parameters that characterize the seeing and AO 
performance. With an estimate of the MTF in 
hand we can compute quantities of photometric 
interest, e.g., the encircled energy versus radius. 

If the distribution of sources in an image is not 
random then it is not a good approximation to 
treat ^(i/)! 2 as a constant and the observed power 
spectrum is not simply proportional to |T(f)| 2 . 
Even if the underlying spatial source distribution 



is random, astronomical luminosity functions are 
commonly steep, and the power in an image can 
be dominated by the the Poisson fluctuations as- 
soicated with a few bright sources. Thus, if we 
do not take account of ^(i/)! 2 when we estimate 
|T(V)| 2 our results will suffer from significant bias. 

We can correct this systematic error in |T(i/)| 2 
by first estimating IS^i/)! 2 from the location and 
relative brightness of stars in the image. This 
can be accomplished by using the methods of 
crowded field photometry to build a model image, 
Smodi*) (e.g., Stetson 1987). The power spectrum 
of s moc ;(x) can be used as an estimate of (S^f)] 2 , 
and the model MTF is multiplied by \S mo d{v)\ 2 
prior to fitting to the observed power spectrum. 
Our results show that even a simple estimate of 
s(x) without any knowledge of the PSF is effec- 
tive. This tactic works well, even when the un- 
corrected seeing halos of multiple adjacent stars 
overlap, which will be a common occurrence in 
AO data. Moreover, this procedure can be applied 
iteratively with DAOPHOT or StarFinder to re- 
fine the PSF and thereby improve s mo( ;(x). If the 
density of stars is sufficient that their diffraction- 
limited cores are blended, then it is unlikely that 
we can measure s m od(x) with any reliability. How- 
ever, this is the circumstance when the approxima- 
tion |5(i/)| 2 = constant may be adequate. 

Phase information has been discarded in Eq. 7 
and some information about the azimuthal struc- 
ture of the PSF has been lost. However, our goal 
here is to estimate the angular scales on which 
is light is scattered so that we can recover the 
photometric curve of growth. Thus, for the cur- 
rent application we choose to consider only the az- 
imuthally averaged PSF and its corresponding 1-d 
power spectrum. Astronomical AO systems do ex- 
hibit non-axisymmetric aberrations — and Eq. (7) 
does encode information about even-order aberra- 
tions, e.g., astigmatism. Thus, modeling the 2-d 
power spectrum may be promising for describing 
effects such as anisoplanatism. 

The data are generally multiplied by a window 
function that represents the field of view and by 
an apodizing function, which is necessary to re- 
duce spectral leakage. In the Fourier domain, the 
window function convolves the Fourier transform. 
Consequently, |/(i^)| 2 is broadened relative to the 
PSF power spectrum by the window function. For 
a square window of width A, the FWHM of the 



corresponding sine 2 function is 0.8859/A Since A 
is typically several hundred pixels, the broadening 
is slight and is ignored in this treatment. 
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Fig. 1.— Calculation of a model MTF for the 
Keck AO system with parameters ro = 65 cm, 
La = 30.0 m, w = 1.5, and a dm = 56 cm: (a) 
Raw atmospheric power spectrum, |$| 2 (dotted), 
the DM filter function, [1 - H'(n)] 2 (dashed), 
and the AO corrected power spectrum, |<&ao | 2 = 
|$| 2 [1 — H'(k)] 2 (solid), (b) Corresponding struc- 
ture functions from panel (a) for the raw (dashed) 
and corrected (solid) atmospheric power spectra. 
(c) The transfer functions (squared). T ( / >AO is cal- 
culated from the AO corrected structure function 
in (b). T pup is the radially averaged approxima- 
tion of the Keck telescope's pupil MTF. T pix is the 
pixel MTF for the narrow camera of NIRC2 with 
a platescale of O'.'Ol pixel -1 , (d) The model sys- 
tem MTF, T sys = T^ AO T pup T pix , plus a constant 
noise floor. The corresponding PSF is shown in 
Figure 2. 



4. Model Keck AO MTF and PSF 

To explore the potential for measuring the MTF 
from astronomical data, we create synthetic im- 
ages similar to those obtained by an actual AO 
system. Ultimately, we will consider data obtained 
using the NIRC2 camera fed by the Keck Obser- 
vatory's AO system operating in LGS mode. We 
thus construct model AO PSFs and MTFs with 
parameters appropriate for this system. 

For a telescope with a circular pupil and a cir- 
cular central obscuration, T pup is a simple analyt- 



ical function of spatial frequency (e.g., Schroeder 
2000). Although dealing with a circular pupil sim- 
plifies the problem, we chose to use the radially 
averaged MTF derived from the approximately 
hexagonal Keck pupil. We first generated a Keck 
pupil image for the NIRC2 camera operating with 
the LARGEHEX pupil stop using the IDL procedure 
NIRC2PUPIL, which is supplied by the Keck Ob- 
servatory 2 . This pupil map includes the central 
obscuration and the spiders. The 2-d autocorrela- 
tion function of this image is T pup . We then aver- 
age this pupil image in equally spaced radial bins 
to construct a radially symmetric pupil transfer 
function. To express T pup in terms of normalized 
spatial frequency, we define a cutoff frequency by 
setting D t ei = 10.99 m, the diameter of the circle 
circumscribing the Keck pupil. The cutoff angular 
frequency of the system, D te i/X, is subsequently 
defined for a telescope with this pupil diameter. 

To construct the AO corrected atmospheric 
MTF, we adopt a modified Kolmogorov atmo- 
spheric phase power spectrum with a finite outer 
scale, 



|$( K )| 2 = 0.0229 r- 5/3 Ll 1/3 (l + « 2 L 2 ) 



-11/6 



(8) 

(cf. Eq. A5) (Strohbehn 1968). Measurements 
suggest that the outer scale of turbulence, L ~ 
10-100 m (Sorrente et al. 2001; Linfield et al. 
2001). The form of Eq. (8) is such that the pres- 
ence of a finite outer scale leaves the power unaf- 
fected, at k> 2tt/L - 

We obtain the AO filtered power spectrum 
|$^, AO | 2 by adopting the influence function appro- 
priate for the 349-actuator Xinetics Inc. mirror 
employed in the Keck system (Wizinowich et al. 
2004). The properties of this mirror are given 
by van Dam et al. (2004) who approximate the 
sinc-likc influence function as the difference of two 
Gaussians, 

h (p) = o — 2 — {exp[-(p/a DM ) 2 /2] - 

A ~ a DM y 

— exp[-2(p/ira DM ) 2 ] } (9) 

where p is the radial coordinate on the DM, and 
odm defines the spatial scale over which a single 
DM actuator affects the wavefront. The Fourier 



2 http://www2. keck, hawaii.edu/optics/lgsao/softwarc 



transform of this circularly symmetric function is 

H(k) = cxp [-2(7r a DM k) 2 ] - 

-cxp[-(7r 2 CT DM K) 2 /2]. (10) 

To mimic the effects of variable correction, we re- 
place the influence function by H'(k) = wH(k). 
A value of w — implies no correction, w = 1 
yields 6 dB, w — 1.5 yields 12 dB, and w = 2 
gives perfect correction at zero frequency. 

Using Eq. (3) for the filtered atmospheric spec- 
trum with an influence function described by 
Eq. (10), and assuming a modified Kolmogorov 
spectrum Eq. (8), we can calculate the atmo- 
spheric MTF by evaluating the integral Eq. (A4) 
for the structure function and substituting the re- 
sult in Eq. (4). The final system MTF is then 
the product of the partially corrected atmospheric 
MTF with the pupil and pixel MTFs (Eq. 5). The 
corresponding PSF for these parameters is then 
found using Eq. (A10). 

As an example, we adopt values of D te i — 
10.99 m, ro = 65 cm, Lq = 30.0 m, <tdm = 56 cm 
(the spacing of the DM actuators projected onto 
the primary mirror), A = 1.65 /im, and w = 1.5. 
Position in the image plane is converted into pixel 
units for the NIRC2 narrow field camera assuming 
a scale of O'.'Ol pixel -1 . We assume that wavefront 
measurement error, A? in Eq. 3, is negligible, and 
set it to zero. The diffraction-limited PSF is com- 
puted by substituting T pup for T sys into Eq. (A10). 

Figure 1 shows the resultant filtered and un- 
filtercd atmospheric phase power spectra, phase 
structure functions, and associated MTFs. The 
filter function, [1 — H'(k)] 2 in Eq. (3), suppresses 
the power spectrum at low frequencies. Figure 1 
shows that the AO system has no effect on spatial 
scales smaller than the actuator spacing. The in- 
ability of the AO system to sense and correct for 
phase errors above the corresponding spatial fre- 
quency along with its imperfect correction within 
the control bandwidth, prevents the final MTF 
from reaching the diffraction limited MTF. The 
resultant narrow MTF produces a PSF broadened 
from the diffraction-limited case with the familiar 
core-halo structure in PSFs delivered by AO sys- 
tems (e.g., Hardy 1998). The diffraction- limited 
PSF (no atmosphere) and the partially corrected 
case are shown in Figure 2. The PSFs shown are 
constructed by interpolating a highly ovcrsamplcd 



PSF onto a grid appropriate for the NIRC2 cam- 
era. Though even a perfect, diffraction limited 
PSF must be sampled onto a pixel array, for the 
purposes of demonstrating change in Strehl ratio, 
the diffraction limited PSF shown has been con- 
structed without taking the pixel MTF into ac- 
count. 




Strehl Ratio = i B.6% 
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Fig. 2. — Model (azimuthally averaged) PSFs for 
Keck-AO/NIRC2 narrow camera (O'.'Ol pixel" 1 ) 
displayed on a logarithmic grey scale. Both PSFs 
have been scaled so that their peak value is 1, 
i.e flux was not conserved. Top: the diffraction- 
limited PSF. Bottom: the partially-corrected PSF 
for ro = 65 cm, Lq = 30 m, o~dm = 56 cm, and 
w = 1.5. The Strehl ratio is 18.6%. 



5. MTF Fitting and Photometry 

5.1. Construction and fitting of simulated 
data 

Using PSFs generated according to the recipe 
of § 4, we generate mock Keck images from which 
we then extract the MTF by fitting our model to 



the resultant image power spectrum. We evalu- 
ate the robustness of our approach by comparing 
the input and recovered parameters and we quan- 
tify photometric errors by evaluating the discrep- 
ancy between the model and reconstructed PSF 
and EE. 

We consider three model PSFs with tq = 50, 65, 
and 100 cm, and common parameters Lq = 30 m, 
w = 1.5, and cfom = 56 cm. Again, we ignore 
measurement noise by setting A? = 0. To explore 
the effects of crowding, we populated images with 
5, 50, 500, 5000, 50,000, and 5xl0 6 stars in a field 
of view of 4'.' x 4'.' . The stars have random positions 
and normally distributed brightnesses. The final 
step of synthesizing a star field involves convolv- 
ing the array of delta functions with the PSF. The 
heights of the delta functions are scaled to approxi- 
mate the distribution of fluxes present in our data, 
and we add normally distributed noise to the im- 
ages to represent detector read noise and Poisson 
sky noise so that the model images have similar 
signal-to- noise ratios as our observations. Some 
example model images are shown in Figure 3. 




Fig. 3.— Model Keck-AO/NIRC2 images demon- 
strating the effects of different values of r$ and 
numbers of stars, N*. All images represent H- 
band data taken with the NIRC2 camera oper- 
ating at the O'.'Ol pixel -1 platescale, and share 
the common parameters Lq = 30 m, w = 1.5, 
o~dm = 56 cm, and A^ = 0. (a) tq = 50 cm, N* 
= 50; (b) r = 65 cm, N* = 50; (c) r = 1 m, N r 
= 50; (d) r = 65 cm, TV* = 500; (e) r = 65 cm, 
iV» = 5000; (f) r = 65 cm, N* = 5xl0 6 . 



Figure 4 shows the measured power spectrum of 



the A* — 50, ro = 65 cm image from Figure 3(b). 
Also plotted is the model MTF (squared) used to 
generate the PSF multiplied by the input source 
distribution power spectrum, |5(^)| 2 as calculated 
from the array of delta functions used to construct 
the image. The model spectrum has been fit to 
the data using only a multiplicative normalization 
factor and an additive constant to represent the 
noise. Inspection of this plot suggests that it is 
reasonable that the MTF should be recoverable 
from astronomical data. 
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Fig. 4. — Measured, model, and best fit power 
spectra for the model image Figure 3(b) (TV* = 50, 
ro = 65 cm). The open diamonds are the power 
spectrum extracted from the image; the solid grey 
line is the known power spectrum computed from 
the input data, consisting of the model MTF mul- 
tiplied by the source distribution power spectrum 
of the original delta function array, and the con- 
stant noise floor; the dotted black line is the best 
fit power spectrum found after two iterations for 
StarFinder to estimate \S mo d{v)\ 2 . 



Since our goal is to establish the unknown MTF 
for real astronomical data we test the ability of 
a non-linear least squares fit to recover the in- 
put parameters from simulated data. We use the 
IDL procedure MPFIT 3 to perform a Levenberg- 
Marquardt least-squares fit of our model MTF to 
Eq. (7) using \S mo d{v)\ 2 as an estimate of the 



source distribution power spectrum. Since we ra- 
dially average the measured power spectrum, we 
are able to directly measure the error in each fre- 
quency bin by calculating the variance of each en- 
semble of pixels. We found that reliable error es- 
timates were necessary to obtain unbiassed fits, 
especially when the AO correction was poor. 

Model/Best — fit Growth Curves 




3 http://cow.physics.wisc.edu/~craigm/idl/ 
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Fig. 5. — Comparison of the input encircled en- 
ergy curves (diamonds) with the recovered encir- 
cled energy curves (solid) from the best fit MTF 
calculated from Eq. (All) for a sample of model 
images, (a) A* = 50, r = 50 cm; (b) A* = 50, 
r = 65 cm; (c) A* = 50, r = 1 m; (d) A* = 5, 
r = 50 cm; (e) A* = 500, r = 50 cm; (f) A* = 
5000, r = 50 cm. 



The model function for T sys (v n ) in Eq. 5 is cal- 
culated in steps. The transfer function T pup (v n ) is 
calculated numerically from the Keck pupil as de- 
scribed in § 4. The pixel transfer function T V i X (v) 
is calculated analytically from Eq. A9. Both T pup 
and T p i X have no free parameters. T^Aoiy) is, 
itself, calculated in steps. We first calculate the 
AO filtered atmospheric phase error power spec- 
trum |$ao( k )| 2 as given in Eq. 3, using Eq. 8 
for the unfiltcrcd atmospheric power spectrum and 
Eq. 10 in H'(n). Calculating these quantities in- 
troduces the fit parameters ro, Lq, &dm, w, and 
Acf,. The structure function is then constructed 
numerically using Eq. A4. (The numerical calcu- 
lation of the structure function is by far the most 
computationally time consuming step, a problem 
compounded by the many iterations necessary of 
a non-linear least squares fit.) T^Ao{ v n) is then 



given by Eq. A7. 

We consider (Jdm & known quantity and hold it 
constant at 56 cm. We also neglect any wavcfront 
measurement error and set A^ = 0. In addition, 
we fit for the multiplicative normalization constant 
of the term \S mo d{v)\ 2 \T sys (v)\ 2 and the noise in 
the image, |A| 2 . These terms, however, do not 
affect the resulting PSF. 

We estimate \S mo d{v)\ 2 with a variety of meth- 
ods. The first uses aperture photometry: for each 
image, we perform fits using the DA0PH0T routine 
FIND implemented in IDL to locate sources in the 
image and to estimate their relative fluxes. The 
second method is iterative: using the parameters 
determined from the FIND fits we generate a PSF 
and use this PSF with the crowded field photom- 
etry package StarFinder (Diolaiti ct al. 2000) to 
again estimate the source function, hopefully more 
accurately. (Note, we use StarFinder only to de- 
termine the location and relative intensity of point 
sources, not to estimate the PSF.) We use param- 
eters from this fit to generate a new PSF to use 
again with StarFinder in a subsequent iteration. 
The third method is used to investigate how well 
the fits work when analyzing model data and the 
source function is known perfectly: we compute 
|5' mo d(^)| 2 directly from the array of delta func- 
tions that was used to construct the model im- 
age. In either case we create an array of size equal 
to the image and place delta functions in this ar- 
ray at the locations of the sources with heights 
equal to the fluxes. We radially binned and aver- 
aged the resultant 2-d power spectrum to obtain 
|<S TOO d(^)| 2 . I n order to leave the overall measured 
power unchanged, we normalize \S mo d{v)\ 2 to have 
a median value of 1. \S mo d{v)\ 2 is then an indi- 
cation of how much ^(j/)! 2 deviates from being a 
multiplicative constant. Finally, we also perform 
fits ignoring the correction for the source distribu- 
tion power spectrum by setting \S mo d{v)\ 2 = 1. 

Figure 4 shows the best fit to the power spec- 
trum for image Figure 3(b) (TV* = 50, ro = 65 cm), 
after two StarFinder iterations. The final fitted 
MTFs are then used to compute the PSFs and 
encircled energy using Eq. A10 and All. Some 
comparisons of the input and recovered EE curves 
are shown in Figure 5. 



5.2. Comparison of input and recovered 
parameters 

In general, the fits are insensitive to the model 
parameters for tq > 150 cm. For tq > 2 m, the 
MTF becomes barely distinguishable from the per- 
fect MTF. In this regime it becomes impossible to 
contrain any of the parameters. It is not surpris- 
ing, therefore, that the fitted parameters for tq = 
100 cm are less accurate than those for r$ — 50 cm 
and 65 cm. This degradation in accuracy applies 
to the fitted parameters for all of the images, as re- 
flected in their calculated formal errors, regardless 
of the number of stars in the field. However, in the 
case of high r , the encircled energy and PSF are 
less dependent on the model parameters, so the 
quantities of interest are left mostly unaffected. 

For images containing 500 stars or fewer, the re- 
covered photometric growth curves deviate from 
the true values by less than 2%. The fits for 
the images containing five and fifty stars gener- 
ally returned the input parameters to within the 
calculated error. The recovered encircled energy 
curves for these fits matched the input well, with 
the maximum deviations from the known growth 
curves less than 1.5%. The accuracy was much 
better than this for ro = 65 cm and 1 to. Estimat- 
ing the source distribution power spectrum and 
accounting for it in the model improved the fits 
and was robust. Using simple DAOPHOT aper- 
ture photometry routine for this purpose was suf- 
ficient, as the separation between the stars was 
large. Signal to noise ratio plays a significant role 
in the accuracy of the fits. The five star images 
yielded less accurate fit parameters than the fifty 
star images. Since we did not hold the total flux 
constant when constructing images with different 
numbers of stars, we created additional five star 
images scaled to match the total signal in the fifty 
star images. Performing the fits on these new im- 
ages yielded more accurate results, in line with 
those obtained from the fifty star images. 

It is more difficult to estimate the source distri- 
bution power spectrum for the images with more 
than fifty stars because of increased crowding. 
Nonetheless, fitting the 500 star images returned 
parameters that were close to the input values 
when we used the iterative method to refine the 
estimate of the source distribution function. The 
resulting photometric growth curves, though less 
accurate than those for five and fifty stars, do not 



deviate from the known ones by more than 2%. 
Using the known source distribution function for 
these synthentic data yields improved fit param- 
eters, confirming that with this many stars un- 
certainty in \S mo d{v)\ 2 is the dominant source of 
error. 

For the images containing > 5000 stars, the ac- 
curacy of the returned parameters drops again. 
Nonetheless, the best-fit photometric curves of 
growth differ from the true ones by usually less 
than 5%. When a field contains so many stars that 
every beam contains multiple stars, estimating the 
spatial distribution power spectrum becomes im- 
possible. In this case, we hope that the assump- 
tion of a random distribution of point sources lead- 
ing to a simplification of ^(V)) 2 as a constant be- 
comes increasingly valid. For the images contain- 
ing 5000 stars (~ 0.4 stars per (A/D) 2 beam), the 
best fit parameters deviate significantly from the 
input parameters. The field is sufficiently crowded 
that errors in \S mo d{v)\ 2 are large, but not so 
crowded that |S'(^)| 2 may be approximated as a 
constant. The deviations in the fitted encircled 
energy from the true encircled energy are as great 
as 5%. These same problems apply to the images 
with 50,000 stars (~ 4 stars per beam), but to a 
greater extent. In this regime, the fits almost com- 
pletely fail to return the input parameters. Accu- 
rately estimating the source distribution here is 
almost impossible. The fits performed while con- 
sidering \S(i>)\ 2 as a constant do not produce bet- 
ter results, indicating that the field is not random 
enough to allow this assumption. 

The 5 x 10 6 image (~ 400 stars per beam) de- 
fies any attempt to estimate the source distribu- 
tion function. However, in this limiting case, the 
field is so crowded that neglecting to account for 
\Smod{v)\ does not cause the fits to fail. The fits, 
though not as accurate as those for the images 
with 5-500 stars, yield encircled energies no more 
than 6% different from the true curves. 

Though poor sky subtraction degrades the ac- 
curacy of the fits, the effect can be mitigated by fil- 
tering out the low spatial frequency bins. We esti- 
mated the background directly from the model im- 
ages and sky subtracted before finding the power 
spectrum. Since the sky subtraction will never 
be perfect and is, for the images with more than 
500 stars, significantly inaccurate, the power spec- 
trum is contaminated. However, if we assume that 



the background does not contain much power be- 
yond the lowest few frequency bins, we can ex- 
clude these bins from our fits. (The first two fre- 
quency bins in Figure 4 demonstrate the effect of 
imperfect sky subtraction.) For the 5, 50, and 
500 star images we ignored the first two frequency 
bins (y. n < 0.02), which yielded more accurate re- 
sults than using all data points. For the 5000 star 
image, we repeated the fits, filtering out the first 
three bins. This resulted in best-fit parameters 
somewhat closer to the inputs, but did not sig- 
nificantly improve the accuracy of the best fit EE 
curve. 

In general, it appears that the benefits of in- 
creased signal to nosic in the images with large 
numbers of stars, as well as the benefits of the in- 
creased validity of the assumption of randomness 
in the source function, are outweighed by the in- 
creased ability to sky subtract and to estimate the 
source function in images with less crowding. (As- 
suming, of course, that there is enough signal to 
be able to use fewer stars.) For instance, the mea- 
sured power spectra for the A» = 5 x 10 6 images 
are noisier than for those with 5, 50, and 500 stars, 
and the fitted parameters are consequently less ac- 
curate, even though the assumption of a random 
spatial distribution holds better for these images 
than for those with less stars and the signal to 
noise is much higher. 

While our most crowded images (A* = 5x 10 6 ) 
yielded fits less accurate than our least crowded 
images (A, < 500), they still led to best-fit curves 
of growth that were accurate to < 6%. There ap- 
pears to be a regime, however, in which the field 
is not crowded enough to allow for the approxima- 
tion of IS 1 ^)! 2 as a constant, but too crowded to 
allow for an accurate estimation of it. Our exam- 
ple of this is the failed fits for images containing 
50,000 stars. In this case, we were unable to find 
a satisfactory way to improve the fits. 

5.3. Fitting laser guide star Keck AO data 

The nearby dwarf starburst galaxy NGC 1569 
(vheiio = —104 km s _1 ) was observed with the 
Keck LGS/AO and the NIRC2 camera on 2005 
March 02. The field is centered on the super- 
star cluster SSC A1/A2 (Origlia et al. 2001). The 
NIRC2 pupil stop was in the LARGEHEX position 
(circumscribed circle) and the camera was oper- 
ating in the narrow field mode (O'.'Ol per pixel). 



Observing conditions were photometric, and the 
seeing was 0'.'4 at H. The R ~ 11 mag. equivalent 
laser guide star was projected at the center of the 
science camera field of view and SSC A1/A2 was 
used as the tip-tilt reference. The science field was 
observed at an airmass of 1.5 with the AO system 
running at 450 Hz. We use a 60-second iJ-band 
exposure from this observation to illustrate the ap- 
plication of our PSF estimation method. Addi- 
tional observations with 120-s exposures were also 
obtained yielding a mosaic with a total exposure 
time of 15 minutes. 
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Fig. 6.— A 60 second Keck LGS AO NIRC2 H- 
band image of super-star cluster A in NGC 1569. 
An MTF was measured and fitted for the entire 
image and for the two indicated subarrays. 



Flat-fields were taken at the zenith near mid- 
night and median filtered to remove starlight. Bad 
pixels in the single 60-s image were repaired by 
substituting them with an interpolation of the sur- 
rounding pixels. As the crowding in the field is 
relatively low and displays little variation across 
the field, conventional sky subtraction was applied 
using the IDL implementation of the DAOPHOT 
routine MMM. 

Figure 6 shows the entire NIRC2 field of view. 
The image is dominated by the resolved cluster 
but also includes several hundred unresolved stel- 
lar sources. The cluster light includes extended 



Fig. 7. — The measured power spectrum of the en- 
tire image in Fig. 6 and the best-fit model power 
spectrum based on the MTF and source distribu- 
tion function, a is the rms error of the ensemble 
of pixels that was radially binned and averaged to 
produce the 1-d spectrum. The fit parameters are 
listed in the last column of Table 1. The corre- 
sponding encircled energy curve is shown in Fig. 



nebular emission and is therefore not suitable for 
measuring the MTF. Our goal is to estimate the 
PSF for the regions containing stars. Therefore, 
we define two subarrays within Fig. 6 that avoid 
SSC A and contain approximately 190 and 20 
point sources, respectively. (Preliminary photom- 
etry indicates that the point sources in the im- 
age have absolute iJ-band magnitudes consistent 
with single red supergiant stars.) Anisoplanatism 
is evident in Fig. 6. We have chosen regions small 
enough (~ 2" x 2" and ~ 4" x 4" ) so that this ef- 
fect should be negligible. For comparison, we also 
compute the power spectrum for the entire image 
after masking out the central cluster and two other 
bright, unresolved objects. 

The best-fit power spectrum and encircled en- 
ergy for the entire image are shown in Figures 7 
and 8. The best-fit parameters are listed in Ta- 
ble 1 along with those for the fits using the two 
subarrays. We again hold a dm at a constant 
56 cm. All parameters represent the best-fit af- 
ter two iterations of using StarFinder to estimate 
the source distribution power spectrum. We also 
state the formal errors in the parameters as com- 
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NGC 1569 Curve of Growth 
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Fig. 8. — The recovered encircled energy curve 
corresponding to the fitted power spectrum (see 
Fig. 7) of the entire image in Fig. 6. The fit pa- 
rameters are listed in the last column of Table 1. 



puted from the covariance matrix. To account for 
any imperfect sky subtraction, we do not fit data 
at v n < 0.02. Thus we filter our spectra below the 
first 1, 2, and 5 frequency bins in the small and 
large subarrays, and in the entire image, respec- 
tively. 
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Fig. 9. — The PSF generated for the entire image's 
best- fit parameters listed in Table 1. 



In § 5.1 our model PSFs were circularly sym- 
metric and based solely on our model MTF. All we 



verified there was the ability to successfully extract 
the MTF from data. Since real PSFs will contain 
asymmetries and since our model will not describe 
all of the factors contributing to an actual MTF 
(i.e., deviations from Kolmogorov turbulence, di- 
rectional dependence in D(r), systematic errors in 
the wavefront sensor) we cannot expect our real 
fits to be as accurate as our fake fits. Nevertheless, 
Figure 7 shows excellent agreement between mea- 
sured and best-fit model power spectra — we have 
evidently included enough detail in our model to 
capture the essence of the Keck AO system. Based 
on the analysis of our synthetic data, the recov- 
ered values put us in the regime where the fits 
show good sensitivty to r$ and w. Additionally, 
the three sets of fit parameters mostly agree to 
within the formal errors. The values of Lq show 
some discrepancy, but in this region of model pa- 
rameter space the MTF and EE show no signf- 
icant variation when Lq changes by as much as 
10%. The recovered EE curve is most sensitive to 
changes in ro, which stays within error over the 
three fits. The density of stars in our image is low 
enough to enable accurate sky subtraction and es- 
timation of the source distribution function. All 
of these factors give us confidence in our derived 
parameters. 

The formal errors decrease significantly as the 
size of the subarray we use to measure a power 
spectrum grows. For one, the signal to noise ratio 
of the power spectrum increases when there are 
more stars. Also, larger images give finer sampling 
of the power spectrum. For these reasons, and 
since the best-fit power spectrum for the whole 
60-s image is a good quantitative match with the 
measured one and does not differ significantly from 
either two subarrays, we choose the entire image 
best-fit MTF to compute a PSF and EE curve. 

We use the best-fit parameters to generate an 
MTF and compute the corresponding PSF using 
Eq. (A10). The recovered PSF is shown in Fig- 
ure 9. This PSF exhibits the same core-halo struc- 
ture that is evident in the data. We performed 
PSF fitting and subtraction for a number of bright 
sources in the image. A few of the original stars 
and the residuals after PSF subtraction are shown 
in Figure 10. The residuals contain approximately 
±10% of the original starlight, a failing that can 
probably be attributed to our initial assumptions 
about the symmetry of the PSF. Indeed, most of 
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Fig. 10. — Typical residuals after performing PSF 
subtraction on stars in Figure 6 using the PSF in 
Figure 9. All images are displayed on a linear scale 
and have a l'/50 x l'/50 field of view. The residu- 
als show mainly azithumal structure, an effect we 
may attribute mainly to the failings of our model 
to describe non-radially symmetric PSFs. Such 
residuals should have little effect on the encircled 
energy. 



the residuals exhibit asymmetric structure. This 
indicates that the actual PSF's deviation from cir- 
cular symmetry, something for which our model is 
unable to account, is responsible for much of the 
residual starlight. 





JV, ~ 20 


N* ~ 190 


Full Image 


r (m) 
L a (to) 

w 


0.58 ± 0.09 
17.5 ± 1.0 
1.33 ± 0.08 


0.56 ± 0.04 
20.0 ± 0.6 
1.34 ± 0.04 


0.56 ± 0.02 
18.1 ± 0.2 
1.33 ± 0.02 



Table 1 : Best fit parameters for the two sub-arrays 
indicated in Figure 6 and for the entire image. 



5.4. Photometry 

In this section we investigate the accuracy of 
photometry extracted from LGS/AO observations 
by comparing the results with data obtained us- 
ing Hubble Space Telescope/NICMOS. The PSF 
for HST instruments is well characterized (e.g., 
Krist et al. 1998) hence the aperture corrections 
are known and HST data provide a good stan- 
dard against which ground-based AO measure- 
ments can be compared. 

Ideally, we would perform the LGS/AO pho- 
tometry in two steps. First, perform relative pho- 
tometry using PSF fitting with an empirical 2-d 
PSF, which is built from stars judged to be iso- 
lated. Second, use our 1-d PSF to compute the 
aperture corrections and place the relative pho- 
tometry on an absolute scale. Since we are inter- 
ested in the accuracy attainable with our method 
we have chosen to perform PSF fitting photome- 
try with the model 1-d PSF obtained as described 
in § 5.3. While crowding will limit the precision of 
this approach, there is no ambiguity in how to cal- 
culate the aperture correction. When applied to 
an ensemble of stars this one-step method should 
give a reliable way of evaluating the accuracy of 
the resultant photometry. 
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Fig. 11.— The Keck AO image of NGC1569A 
and the corresponding sub-image of the NICMOS 
mosaic. The NICMOS image has been rotated to 
match the Keck field-of-view. The arrow indicates 
a variable star with the largest observed fluctua- 
tion (5m = 0.8 mag.). 



We measure iJ-band magnitudes for approxi- 
mately 250 bright point sources in the field of a 
single 60-s exposure of NGC 1569 SSC A (r = 
56 cm; Strehl = 10.1%; see § 5.3). We used five 
images of the near- infrared standard star SJ 9139 
(Persson et al. 1998) obtained on the same night 
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as the science data to measure the photometric 
zero point. 
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Fig. 12. — The magnitudes for 243 stars common 
to the Keck and HST fields of view. The HST 
magnitudes were taken from Aloisi et al. (2001) 
and color corrected from the NICMOS F160W fil- 
ter to Keek's H filter. The dashed line has slope 1 
and intercept 0. The triangle is the star indicated 
by the arrows in Figure 11. 



appendix of Marleau et al. (2000). 

A plot of Keck vs. color-corrected HST mag- 
nitudes for isolated stars (Figure 12) shows that 
the mean photometric discrepancy between the 
two data sets amounts to only a few percent. In 
this figure the data points fall along a straight 
line of unit slope and intercept zero. The mean 
Keck magnitude minus the HST magnitude, 5m, 
is 5m — 0.05 mag., and standard deviation is as m 
= 0.21 mag. The standard error of 5m is 0.013 
mag. 
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HST data were obtained with the NIC2 cam- 
era in the F160W filter on 1998 February 25. The 
Keck image and the matching NICMOS sub-image 
are shown side- by-side in Figure 11. Photome- 
try for the HST/NICMOS images has been ex- 
tracted by Aloisi et al. (2001). Comparison of 
the photometry obtained from our LGS/AO data 
begins with matching the catalogs of Keck and 
HST/NICMOS sources. Inspection of Figure 11 
shows that Keck frequently resolves stars that ap- 
pear as single sources to HST/NICMOS. Although 
we are able to identify more stars in the Keck 
image than presented here, we have selected only 
those stars that are separated from neighbors by 
more than O'.'l in both catalogs. The central wave- 
lengths of the NICMOS F160W filter (1.60 /xm) 
and the NIRC2 iJ-band filter (1.63 /zm) are sim- 
ilar. However, the filter bandpasses differ signif- 
icantly, and for these red stars a color correction 
is necessary to compare space and ground-based 
measurements. We correct the HST F160W mag- 
nitudes to the Keck system based on their F160W- 
F110W colors, using the method outlined in the 



Fig. 13. — The difference between the measured 
60-s Keck exposure magnitudes and the HST mag- 
nitudes as a function of distance from the center of 
the FOV. Anisoplanatism would appear as an er- 
ror that increases with distance from the laser bea- 
con (projected at r = 0). The dashed line shows 
the best-fit second order polynomial, illustrating 
that there is only a weak trend in the photometric 



The agreement between the Keck and HST 
photometry at the 5% level is encouraging, and 
provides significant validation for our method of 
finding AO aperture corrections. However, the 
roughly 20% scatter for individual measurements, 
which is several times larger than that expected 
from source, sky and detector noise, merits some 
additional comment. 

One possible source of scatter is anisoplanatism 
in the AO images. A systematic change in the 
PSF and across the field of view would result in 
seemingly random errors distributed among stars 
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of all magnitudes in Figure 12. However, a plot of 
dm vs. field angle (Figure 13) shows only a weak 
trend, and indicates that variation of the AO PSF 
across the Keck field of view only contributes a 
fraction (~ 10%) of the 0.21 mag. scatter. 

We investigate our internal errors by compar- 
ing two independent sets of LGS/AO observations 
of the field of NCG 1569 SSCA. For the first im- 
age we use the same 60-s exposure as above and a 
120-s exposure for the second (tq = 42 cm; Strehl 
= 3.7%. Each image has a different PSF and a 
sky background level. Comparing the photome- 
try done from the 60s and 120s exposures gives 
m 60s - toi2o s = 0.04 mag. and 0- (m6Os _ mi2Os) = 
0.10 mag. Thus, our measurement errors account 
for only about 20% of the variance between the 
Keck and HST. 

Likely, the major source of the observed scat- 
ter is the natural variability of red supcrgiants 
(RSGs), the main stellar component present in 
this field. Visual inspection of the NIRC2 and 
NICMOS images reveal several stars that have 
brightened or dimmed betweeen 1998 and 2005. 
The most variable star, with 5m = 0.8, is indi- 
cated by the triangle in Figure 12 and the arrow 
in Figure 11. The magnitudes and F160W-F110W 
colors (« 1) of our sample indicate they are mostly 
RSGs with masses > 10M Q (Aloisi et al. 2001). 
The most massive few stars are probably between 
20M© and 3OM0 . The stars in our sample arc un- 
stable and can fluctuate significantly in luminos- 
ity. Maeder (1980) has studied this variability in 
a sample of 2420 RSGs. The results indicate that 
stars of these masses and spectral types have av- 
erage peak-to-trough U-band amplitudes between 
0.1 for 10M Q stars and 0.4-0.8 for 3OM stars. Al- 
though the amplitude of variation is likely to be 
less at infrared wavelengths, these findings suggest 
that we can attribute most of the observed vari- 
ance to the stars themselves. 

Since the 120-s exposure exhibited poorer im- 
age quality, it affords the opportunity to examine 
the photometric integrity for reduced AO correc- 
tion. The same photometric Keck-HST compar- 
ison for these data give 8m = —0.03 mag. and 
c<5m = 0.23 mag., confirming that accurate results 
can be extracted at low Strehl ratios. 



6. Discussion 

We have developed a simple model to describe 
the imaging performance of an AO system. The 
model makes assumptions about the shape of the 
power spectrum of atmospheric phase errors, the 
pupil geometry, and the DM influence function 
and has a few adjustable parameters, including 
the Fried parameter (ro), the degree of AO cor- 
rection (w) , and the wavefront measurement noise 
(Ac/>). Beyond the inherent value of a well defined 
framework to quantify an observed phenomenon, 
this method has several advantages in the practical 
analysis of AO data. Generating numerical MTFs 
is a cheap and easy way to predict system perfor- 
mance under a variety of conditions, i.e., seeing, 
pupil obscuration, DM actuator spacing, &c. 

By considering the Fourier information present 
in astronomical images we have shown that it is 
practical to fit model MTFs to image power spec- 
tra, constrain the model parameters and recover 
the PSF, encircled energy and Strehl ratio. High 
fidelity simulations of AO performance can be 
performed using computationally intensive Monte 
Carlo codes (e.g., Jolissant 2004). However, we 
believe that we have shown that it is unnecessary 
to employ these elaborate methods to compare 
with observations to estimate photometric curve 
of growth corrections. 

Using the method outlined here, the az- 
imuthally averaged PSF and hence the photo- 
metric curve of growth can be estimated directly 
from an image containing multiple point sources. 
Data obtained in varying observing conditions 
and with changing AO correction can be ana- 
lyzed so that they are not subject to systematic 
photometric errors, and separate and potentially 
unreliable PSF star observations are not required. 
Our method allows one to estimate the PSF, even 
if there are no isolated stars in the image. Al- 
though only the azimuthally averaged PSF is de- 
rived, the fact that the residuals in Figure 10 
show mainly azimuthal components means that 
the curve of growth is valid. A comparison of 
Keck LGS/AO and HST/NICMOS observations 
of the same field at 1.6 ^m in NCG 1569 shows 
that our photometric method is accurate to 5% 
or better. The precision of our 1-d PSF fitting 
photometry (about 10%) is not limited by photon 
statistics and detector noise. We are investigating 
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a hybrid approach which involves conventional 2-d 
PSF fitting photometry using an empirical kernel 
followed by application of aperture corrections de- 
rived from MTF fitting. 

Treating the OTF as a one-dimensional MTF is 
an approximation. The wavefront delivered by an 
AO system often suffers from systematic calibra- 
tion errors which lead to significant PSF azimuthal 
asymmetry. Non-circular pupils and anisotropic 
atmospheric turbulence, e.g., due to a prevail- 
ing wind, are the rule rather than the exception. 
While odd-order aberrations, e.g., coma, cannot 
be reconstructed from the power spectrum, even 
aberrations like astigmatism can be described. A 
natural extension to the work presented here is to 
fit a 2-d OTF to the 2-d power spectrum to cap- 
ture this information. It will be particularly in- 
teresting to investigate the ability of this method 
to describe the astigmatism introduced by angu- 
lar anisoplanatism by mapping out the variation 
of model parameters as a function of field angle. 
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A. Appendix: Modulation Transfer Functions 

The optical transfer function (OTF) specifies the amplitude and phase response of a linear optical imaging 
system to a stimulus at a given spatial frequency. In general the OTF is complex. Here we consider the 
special case of a circularly symmetric PSF, in which case the complex OTF reduces to the real modulation 
transfer function (MTF). Following Schrocdcr (2000), the MTF, T(v), is a measure of the contrast at spatial 
frequency v in the image, relative to the object. We expect T(y) — > 1 as v — > 0. For a pupil with finite 
support, above some critical spatial frequency, v c , T(v c ) — and all information at frequencies greater than 
v c is lost. 

In the image plane, the cutoff spatial frequency of a telescope with a circular pupil, diameter, D te i , and 
observing wavelength, A, is 

v c = 1/FX. (Al) 

This corresponds to a linear image size of f\/D te i — FX, where / is the effective focal length at the image 
plane and F = f / D te i is the focal ratio. It is convenient to use spatial frequencies in units of the cutoff 
frequency in the image plane, or normalized spatial frequencies, v n — v jv c = v\F . 

A.l. Atmospheric MTF 

The variation in phase of a wavefront, (f>, between two points in the pupil plane of a telescope separated 
by a distance r = ri — r 2 can be described by the structure function 

^ (r)^([0( ri )-0(r 2 )] 2 ), (A2) 

where (...) denotes a time average. We reserve lower case, e.g., </>, to refer to quantities in the spatial domain, 
and upper case, e.g., <J?, for the Fourier transformed quantity in the frequency domain. We use x and v 
to refer to spatial coordinates and spatial frequencies in the image plane and r and k to refer to spatial 
coordinates and spatial frequencies in the pupil plane. Thus, the relation between position in the pupil plane 
and spatial frequency in the image plane is r = Xfv. 

The structure function can be written in terms of the covariance function, which in turn forms a Fourier 
pair with the two-dimensional phase-error power spectrum, |$(/t)| 2 , where k is spatial frequency in the pupil 
plane. The notation |<I>| 2 denotes the atmospheric phase error power spectrum, which is shorthand for the 
expectation value ($$*), where * denotes the complex conjugate. Consequently, the structure function and 
the power spectrum are related by the transform 

D^r) = 2 I |$(/t)| 2 [l - cos(2ttk • r)]d/e, (A3) 

where the integral is taken over the infinite domain. In the case of an isotropic power spectrum integration 
over the angular part yields 

c2tt 



/>00 />Z7T 

D (r) = 2/ \$(k)\ 2 Kdn d9 [1 - cos(2ir kr cos 0)} 

Jo Jo 

/>oo 

= 4tt / |$(k)| 2 [1 - J (2ir Kr)]K dn 7 (A4) 



'o Jo 

JO 

where k = \k\ (Tatarski 1961), and Jo is the Bessel function of the first kind. Note that this transform, as 
quoted by (Hardy 1998, Eq. 4.46), contains a typographical error. 

For Kolmogorov turbulence, with an infinite outer scale and an inner scale of zero, the two dimensional 
phase error power spectrum in the pupil plane is 

|$(k)| 2 = 0.0229 r(7 5/3 K- n / 3 (A5) 

(Noll 1976). The corresponding phase structure function is 
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/ \ 5/3 
D^(r)= 6.88 ( -J , (A6) 

where r is the separation between two points in the pupil plane and tq is a scale length known as the Fried 
parameter (Fried 1965). The phase structure function evaluated at r — Xfv — D te iv n yields the atmospheric 
MTF structure function 

T (i/„) = cxp [-\D^D tel u n )] (A7) 

(Fried 1966). 

A. 2. MTF of an unaberrated pupil 

In the case of a unaberrated, unobscured, circular pupil, the OTF reduces to the real, circularly symmetric 
MTF, 

T pup {v n ) = - [cos" 1 v n - i/ n (l - v 2 S l2 \ . (A8) 

IT L J 

In the more general case of a pupil function with central obscuration, secondary spiders and segment gaps, 
the MTF can be computed numerically from the autocorrelation of the pupil function. 

A.3. MTF of an ideal pixel 

The MTF of a ideal pixel of dimensions a x a can be derived by considering it as a uniformly illuminated 
square aperture. The MTF is then given by the Fourier transform of this function, which in normalized 
frequencies is 

Tpi x (v n ) = sinc(7n/„/n p ix), (A9) 

where n p j x is the number of pixels per length XF. For Nyquist sampling, two pixels per XF, n p ; x = 2. For 
the narrow field camera in Kcck/NIRC2 n p ; x = 3.4 at H band. 

A. 4. PSF, encircled energy, and Strehl ratio 

The final step involves finding the PSF, p(r n ), and the encircled energy, e(r n ), from the Hankel transforms 
of the system MTF T ays (i/ n ) = T^(fc'„)Tj, U p(z/„)Tpj a: (i/„). The Hankel transform of zero order yields the PSF, 

p(r n ) = 2n T sys (v n )Jt)(2iTv n r n )v n dv n (A10) 

Jo 

where r n = r/FX = 9D te i/X is the normalized distance in the image plane, and the Hankel transform of first 
order yields the the encircled energy 

e(r„) = / p(r' n ) 2nr n dr' n = 2nr n / T sys (v n ) Ji(2irr n v n ) dv n . (All) 

Jo Jo 

where Jo and J\ are the Bessel functions of the first and second kind, respectively. The Strehl ratio is defined 
as Psys(0)/p P u P (0), where p pup denotes the PSF for the unaberrated system. Thus, using Eq. (A10) 

»i 
T sys (u n )v n du n 
S=^- . (A12) 

J- pup \ l^n ) ^n wl^n 
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